Absence of Dipole Transitions in Vortices of Type II Superconductors 
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The response of a single vortex to a time dependent field is examined microscopi- 
cally and an equation of motion for vortex motion at non-zero frequencies is derived. 
Of interest are frequencies near / Ep, where A is the bulk energy gap and Ep is 
the fermi energy. The low temperature, clean, extreme type II limit and maintaining 
of equilibrium with the lattice are assumed. A simplification occurs for large planar 
mass anisotropy. Thus the results may be pertinent to materials such as NbSe2 
and high temperature superconductors. The expected dipole transition between core 
states is hidden because of the self consistent nature of the vortex potential. Instead 
the vortex itself moves and has a resonance at the frequency of the transition. 
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The description of quasiparticle levels in vortex cores of superconductors has been known 
for some time. Caroli et al. and Bardeen et al. calculated the energies and wavefunc- 
tions of these discrete levels using the Bogoliubov-deGennes (BdG) equation. This was a 
basis for theories of dissipative vortex motion based on the idea of a 'normal' core such as 
those of Bardeen and Stephen 0, Nozieres and Vinen (NV) and others Kramer and 
Pesch p used the Eilenberger equation to calculate the density of states in the core. Their 
approach has proven useful in qualitatively explaining [0 scanning tunneling microscope 
experiments on NbSe2 by Hess et al. These experiments, however, probe the density of 
states at a scale ~ O.lmeV whereas the separation of core levels is ~ lOmK. Caroli and Ma- 
tricon discussed the implications of discrete levels for ultrasound attenuation and nuclear 
magnetic relaxation. Transitions between these levels may have been observed recently in 
high temperature superconductors |T^. The electromagnetic response is interesting from a 
practical point of view. Herein we focus on the low temperature and clean limit, considering 
eigenfunctions and matrix elements of quasiparticle states, and assuming the BdG equations 
and a local gap equation, A (r) = V{ci (r) C| (r)), {c^ is a spin a electron operator) to be 
valid. The vortex response to an electromagnetic field will be considered from a purely 
microscopic point of view. Real materials have vortex pinning but at frequencies of order 
of magnitude comparable to the core level separation undergo a cross-over to unpinned be- 
haviour so a first step should be a study of unpinned vortices at those frequencies. Real 
superconductors are non-local and the interaction is retarded but the hope is that, since 
we deal with rigid motions and relatively low frequencies, some relevance to real materials 



remains. 



In terms of eigenfunctions, ilj^{r) = {u^{r) v^{r)), the quasiparticle operators are, 
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The Schrodinger equation for ip is the BdG equation. 
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A*(f) 

where A = |A(r — fo)| exp {—i6{f— fo))- ^(^— ^o) is the angle about the center of the vortex 
To measured from the x axis. We consider a vortex parallel to z. 

In the extreme type II limit with H « Hc2 the magnetic field creating vortices may 
be ignored. Its importance compared to the phase of A is reduced by where ^ is the 

coherence length and A is the penetration depth. 

The eigenfunctions for fixed k^, fi << kp±C,, and the radial coordinate r << ^ are 
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where /i = ±i,±|,... and ± refers to the x,y directions (we assume at least cylindrical 
symmetry). They fall off exponentially for r > C,- The energies as calculated by Kramer and 
Pesch P], who accounted for some self-consistency effects due to the gap equation, are 



= , ^^^°o^ ^^(^^ocog6/^i), cosQ = kF±/kF. 



(0.4) 



The logarithmic factor is not important here and ignored hereafter but see reference [12 



Consider a long wavelength electromagnetic wave, A ± 5, with polarization at angle 
to X. We treat the perturbation — —A ■ p to first order in A'. If the matrix elements with 
respect to low energy states are not very sensitive to the exact behaviour of the wavefunctions 
at the core boundary then they may be estimated, using standard Bessel function identities, 

to be ■V)^^^ = {ehkF±A'/2mc)expT^{0o + f)- 

Consider a vortex with velocity vl -i- z a. background superfluid velocity vs- These 
velocities will be assumed uniform along the length of the vortex. This is valid if the 
distance the electromagnetic wave penetrates the superconductor (the shorter of the London 
penetration depth and skin depth) is long compared to the coherence length. We shall 
consider a background supercurrent or gauge field of the form. 



A'{t) = -{mc/e)vs{t) = Ect + A'q. 
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E is the applied electric field. 

Let us now outline the ensuing calculation. The time dependence of the quasiparticle 
states under the applied field and the moving vortex will be calculated and written in terms 
of a density matrix. The motion of the vortex will be inferred, using the gap equation, by 
identifying changes in the off-diagonal components of the density matrix with that due to 
a displaced vortex. Because this calculation involves self-consistency in the vortex velocity 
there is a slight subtlety. Given inital (t=0) values of vl and vs, as considered below, in 
general they will not be consistent. In general, vl is obtained by integrating the equation 
of motion for a given vs{t). We shall first derive the acceleration of the vortex at t = 
and then find the equation of motion for all time by beginning in a well defined equilibrium 
state { Vl = 0, Vs = 0, fo = at t = 0) and, given a time dependence vs{t), calculating the 
resulting motion to all orders in t. 



Suppose the field in equation is turned on at t = 0. Let us calculate the quasiparticle 
density matrix to order t"^. We use ()o to denote the (diagonal) values at t = 0. Using matrix 
elements of A' the off-diagonal density matrix elements are, 

{lll,-i)it) = [{lll,)o - {ll-il,-iU X - ^ [ihW + W{e,., - e^)]| (0.6) 

where W = {ehkF±AQ/2mc) expi(6'o -l- f ) and W is the same as W except with A'q replaced 
by cE. Terms of order W"^ (such as the change in the diagonal density matrix element) are 
ignored by taking the amplitude of the perturbation to be sufficiently small. 

The significance of these density matrix elements is clarified by considering a displaced 



vortex in terms of the underlying quasiparticles. The inverse of equation ^3 substituted 
into the gap equation is 

= Vj2[i^,. - E(7L7m-))<(^1«/.(0 + other terms]. (0.7) 

Suppose the quasiparticles are displaced by 6ro at angle 0o to x. The occupation is taken 
to be diagonal before displacement, {"il^l ^j.a) = Sfiufi^u), where /(e) is the fermi function. 
Then, in the undisplaced eigen-basis, to linear order in (5ro, 



SA{r) = -VSrokp^Yl [e"^°ifi^^) ' f ie,+i)H^i{^u,{^ (0.8) 

V 

+ e-*^"(/(e.) - /(e.+i))<(f)«,+i(r-0]. (0.9) 



Comparing with equation p.6| it is clear that the term of order t may be produced by a 
rigid translation of the vortex core at velocity 5ro/t = — (e/mc)AQ in direction 0o = ^o- It 
is simply the velocity by which the gauge field boosts the group velocity of all waves. 

Returning to equation p.6| the first term of order corresponds to an acceleration propor- 
tional to the electric field and equation ^]9| shows that it is vs- The second piece is trickier. 
Because e^_i — can depend on kz and (for large /i) /x, the term does not correspond 
to a rigid acceleration of the vortex core except in the limit of low temperature and strong 
mass anisotropy, >> m_|_. There e^_i — = 2A^/kpVF (independent of yu) and making 
the identification with a rigid displacement leads to an acceleration 



25ro _ 2A2eA[, 
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t"^ hkpVpmc 

in a direction </>o = + f perpendicular to vs- This corresponds to the Lorentz force. 

At high temperatures higher energy levels become important. The level spacing decreases 
with increasing energy so higher energy quasiparticles have a smaller acceleration and lag 
behind the core. The term e^_i — contains a factor (1/cos^G) which is the kz dependence. 
This factor is not actually divergent. The expression |0]J for the energies is valid only for 
small energies. The are bounded by Aq and so e^_i — must go to zero as 6 7r/2. 
In real materials such as NbSe2 the fermi surface is open which restricts cosQ. Below, for 
simplicity, we assume that the system has an anisotropic mass and that the lack of rigid 
acceleration is not important. 

Assuming rigid motion we see that the applied field, instead of causing dipole transitions, 
causes the density matrix to evolve off-diagonal elements corresponding to vortex motion 
(after applying the gap equation). In the new displaced set of basis functions the density 
matrix is again diagonal. The vortex does not stand still and allow a dipole transition to 
take place, as the core of an atom. The vortex is a self consistent potential. 
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The moving vortex itself affects the density matrix. Suppose the vortex has velocity vl 
at an angle (p^io x. The matrix element 



( (5A(t) ^ 



^.(r) (0.11) 



y<5A*(t) J 

may be re- written using Galilean invariance. Let 5il)^{r) be the change in ipnif) upon 
displacement of the vortex by 6fo. Then substituting 6ip and 6 A into the BdG equation, 
subtracting the undisplaced piece, keeping terms to first order in 6rQ and integrating by 
parts results in W^^j = {ty — e^) / ip'^^ij^dipyir). Writing 5ip{f) = tv^ ■ Vip{r), and using 
standard identities, 

= -^^5,,^i(6. - 6,)e±^(^o-^f ). (0.12) 

Treating this to linear order and integrating from to t gives an acceleration (25ro/t^) = 
(yL/ft-){^u — ^u+i) = —{Al/Eph)vi^ in a direction at an angle +n/2 to v^. 

This together with equation p.lO| gives an acceleration {Al/hEp){vL — vs) x z corre- 



sponding to the Magnus force given in reference |Q as, {hn/2){vs — vl) x 5, where n is 
the (superfiuid) electron density. Taking n = kl^^kpz/Ti^ and the in-plane coherence length 
= hvF±/'n'A, one may extract a 'mass' of the vortex M ~ m±{kF±C,±y per unit length 
kpl- This expression is perhaps a microscopic justification for a 'normal core' of size ^, even 
at low temperatures when there is a gap in the single particle density of states. This mass 
should be contrasted with the very different definition of mass discussed recently by Duan 
and Leggett [|1^. Here the mass corresponds to the inertia of the electrons in the core of 
the vortex. 

To treat dissipation we assume that core states maintain equilibrium with the lattice and 
that the linear in t change in the single particle states (in the lattice frame of reference) also 
decays as an exponential exp (— t/r). r is related to the transport lifetime although there 
are differences in matrix element and phase space which can be elucidated in a microscopic 
theory. Given the velocity vl, then to linear order in t the off-diagonal component of the 
density matrix is, in a basis fixed to the lattice at t = 0, 
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= (7^^)(1 - {lUl,-i))[-{kF±VL/2)e^'^'^]t (0.13) 
+ (7i-i7M-i)(l - {lil,)mFLVL/2)e^'^^^]t (0.14) 

The result is an additional contribution to the second derivative of the off-diagonal element, 
(f Id 

^(7i7/.-i) = --^(7i7/.-i) + (previous terms). (0.15) 

This produces a contribution —{l/r)vL to the vortex acceleration. Collecting terms the 
equation of motion is 

1 

vl = vs + T-^-i^L - Vs) X Z Vl- (0.16) 

In the above derivation we applied ^5 instantaneously to a vortex stationary for t < 0. 
That produced vl = vs B:t t = which is an insufficiently general initial condition. We now 
show that equation p.l6| is valid at all times. We begin at t = in a well defined state, 
Vl = vs = 0, and displacement ro = 0. Given the Taylor expansion of vs we may calculate, 
using the above technique, the displacement to any order t" because we need know only vs 
to order n — 1 and vl to order n — 2. The gap equation applied to the order n displacement 
of the quasiparticles gives vl to order n — 1 and thus continues the calculation. The result 
to fourth order is, 

ro(t) = \^sm^ + It (^^5(0) - r-^'vs{^)) (0.17) 
+ i [-T-'nolMO) X ^] - r-%{0) + T"%{0)) + • • • (0.18) 

where Qq = A'^/hEp. This calculation is fully consistent and equivalent to taking the 
derivatives of equation p.l6| and evaluating them at t = 0. Since the equation is linear and 



velocities can be calculated to all orders in t this equation is valid for all t > 0. 

Equation p.l6| at zero frequency was introduced by deGennes and Matricon |T^. The 



dissipation acts on vl rather than vs as in the NV equation. This has the drawback of not 
allowing for small conductivities observed in experiments on flux flow [|T^ (not to mention 
that there has never been a satisfactory explanation of the Hall effect). The present deriva- 
tion is valid in the clean, low temperature limit where the levels are clearly separated. Our 



result agrees with the NV equation of motion in that hmit [jT6[ but cannot be extended to 
the dirty hmit. 

It is easily verified that the homogeneous solution of |0.16| corresponds to circular motion, 
with a definite handedness, at frequency Qq and decaying on a timescale r. To obtain a pre- 
diction for the surface impedance consider the particular solution for vs{t) = Vsi^) exp (iut). 
It is 



VLy - Vsy 



^O'TVsx - {l + iuJT)vsy 



(0.19) 



(l + ^cur)2 + (^]or)2 

and another equation with x,y interchanged and Qq —^o- 

There are two contributions to the surface impedance. The first is transverse vortex 
motion in phase with the supercurrent. For clarity let vsy = Vsx^^pid = Vs/V^- There is 
an induced voltage per vortex z x VL[h/2e). The supercurrent density vsne gives dissipation 
N^{hn/2) Re{v1^v sy — vlyVsx), where is the vortex density. Using equation p.l9| 

ReivLvsy - vlyvsx) = -vl^or (0.20) 

The first term has an (anti)resonance at (u^r)^ = {^lory + 1 while the second, polarization 
dependent, term has a resonance there. 

The second source is the current due to vortex motion which is in phase and parallel with 
the applied electric field. A straightforward calculation gives the average current density due 
to vortex motion to be 2vi{kFzl'^){,EF/ ^^N^e. With the electric field E = —{m/e)vs the 
dissipation is —2m{kFz/T!'){Ep/AyNyRe{vlyiu!Vsy + v*j^J.ujv sx) ■ Using equation p.l9| , 

ReW^ru^vsy + v,^^uvsx) = ^ _ + (^^^pp + 4(^^)2 • (0-22) 

This has a peak at {ujtY' ~ (noT)^ + l and is not polarization dependent. By letting VLq — »■ 0, 
this term becomes simply the Drude expression for dissipation. 

The expected polarization dependent absorption is distributed, due to states with dif- 
ferent A;^, over a range of frequencies. Precise details depend on the fermi surface shape. 
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The temperature dependence should be weak. The current due to quasiparticles moving in 
and out of the vortex core can be neglected at low temperature because there is a discrete 
energy cost to make charge fluctuations in the core. The considerations of this paper may 
be valid even for pinned vortices. If some parts of a line are pinned, other parts of between 
pins can move, provided they are excited at high enough frequency. 
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